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ABSTRACT 

KOI- 13 was presented by the Kepler team as a candidate for having a giant planet — KOI-13.01, 
with orbital period of 1.7 d and transit depth of ~0.8%. We have analyzed the Kepler Q2 data 
of KOI- 13, which was publicly available at the time of the submission of this paper, and derived 
the amplitudes of the beaming, ellipsoidal and reflection modulations — 8.6 + 1.1, 66.8 ± 1.6 
and 72.0 +1.5 ppm (parts per million), respectively. After the paper was submitted, Q3 data were 
released, so we repeated the analysis with the newly available light curve. The results of the two 
quarters were quite similar. From the amplitude of the beaming modulation we derived a mass 
of 10 ± 2M Jup for the secondary, suggesting that KOI-13.01 was a massive planet, with one of 
the largest known radii. We also found in the data a periodicity of unknown origin with a period 
of 1.0595 d and a peak-to-peak modulation of ~60 ppm. The light curve of Q3 revealed a few 
additional small-amplitude periodicities with similar frequencies. 

It seemed as if the secondary occultation of KOI- 13 was slightly deeper than the reflection peak- 
to-peak modulation by 16.8 ± 4.5 ppm. If real, this small difference was a measure of the thermal 
emission from the night side of KOI-13.01. We estimated the effective temperature to be 2600 ± 
150 K, using a simplistic black-body emissivity approximation. We then derived the planetary 
geometrical and Bond albedos as a function of the day-side temperature. Our analysis suggested 
that the Bond albedo of KOI-13.01 might be substantially larger than the geometrical albedo. 

Key words. Methods: data analysis — planetary systems — stars: individual: Kepler KOI-13 — 
KIC 9941662 



1. Introduction 

The Kepler team ( Borucki et ailboi 1 ) listed KOI-13 (KIC 9941662) as a candidate for h aving a gi- 
ant planet with an orbital period of 1 .7 d and a transit depth of ~0.8% dRowe et al.1 (1201 lb discussed 
this system in an abstract). The planet got a Kepler vetting score of '2', which, according to the 
classification of Borucki et al., meant that the planet was a "strong probability candidate, cleanly 
passes tests that were applied". This score was only next to that of the planets confirmed with 
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radial-velocity (RV) follow-up observations, which got a vetting score of ' 1', and w ere then chris- 



tened an official planet sequential Kepler numbering, like Kepler- la, Kepler-2a, etc dBorucki et al 



2010). 



The Kepler Inpu t Catalog (IBrown et al.ll201 II) stellar temperature of KOI-13, as appeared in 
(1201 lb . is 8448 K. Stars with such temperature have very f ew, usually wide, ste llar 



Borucki et al 



Galland et al. 



2005). It 



lines in their spectra and therefore measuring their RV is difficult (e.g., 
seemed therefore that the planet around KOI-13 — KOI-13.01, was doomed to stay with the '2' 
vetting score, and might not rea ch the official Kepler -planet exclusive club. 



This is unfortunate, because 



Borucki et al 



(201 1), based on the KIC stellar mass and radius and 
the analysis of the Kepler light curve, estimated the planet candidate radius to be 1.86 + 0.003/? j up . 
This is one of the largest derived radii of transiting planets. According to the Extrasolar P lanet 



2011) and 



Encyclopedia (http://exoplanet.eu/), as of July 2011 only WASP-17b (lAnderson et al. 
HAT-P-32b dHartman et alj|201 ll) have larger radii, of about 2Rj up . These two confirmed planets 
have masses of 0.5 and 0.9 Mj up , respectively, and therefore are some of the most bloated, low- 
density planets known. It would be interesting to have an estimate of the mass of KOI-13.01 and 
see if it is also one of th ose bloated p l anets. 



In fact, very recently 



Szabo et al 



(201 1) used the derived large radius to argue that KOI-13.01 
was not a planet at all. Szabo et al. noticed that KOTO was the known BD+46 2629 visual double 
star, and convincingly proved, by analyzing the publicly available Kepler images, and by obtaining 
ground-based lucky imaging photometry, that the planet candidate orbited the brighter of the two 
stars — KOTO A. They also showed that the light coming from KOT13B contributed 45% of 
the light of the system in the Kepler band. Furthermore, Szabo et al. obtained a high-resolution 
spectrum of KOT 1 3 and derived its mass, radius, temperature and stellar rotation — v sin i rot , where 
i ro , was the inclination angle of the stellar rotation axis. The latter was estimated to be 65 + 10 
km/s. They also discovered that the transit had an asymmetric part, on the order of 100 ppm (parts 
per million), and suggested that this was caused by some obliquity of the orbit and high gravity 
darkening of the stellar surface, due to its high te mperature and fast ro tation. After correcting for 

(ho 111) derived the radius of KOI- 



Szaboetal 



the dilution factor and the gravity darkening effect, 
13.01 to be 2.2 + 0.1 R j up . They used this value t o argue that KOI-13 .01 was a brown dwarf, as the 
theoretical models of planet evolution they cited ( Fortnev et alj|2007 ) could not account for such a 
large radius. 

The conclusion of Szabo et al. seems somewhat premature, given the fact that we know of HAT- 
P-33b and WASP-17b, two bona fide planets with RV mass determination that have radii ~2Rj up . 
We would guess that the theoretical border of the allowed range of planetary radii is not so sharp. 
It is therefore even more desirable to have some observational estimate of KOI-13.01 mass, so we 
can confront the theory with the observations. 

One way to esti mate the mass of KO I -13.01 is to use the Doppler beamin g (sometimes called 



2003 



Zucker. Mazeh & Alexander 



2007) and the ellipsoidal 



Doppler boosting) (Loeb & Gaudi : 

(see a review by Mazehl20*08 ) effects. Both effects cause a modulation of the stellar brightness with 
the orbital period (beaming) and its first overtone (ellipsoidal). To first approximation, both effects 
depend linearly on the secondary mass and therefore detecting these two effects and measuring 
their amplitudes can, in principle, yield a ma ss estimate for the planet candidate. Such an approach 
was applied recently bv ICarter et al.l (1201 ll) to KIC 10657664 — an A-star with a white-dwarf 
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secondary, although the results from the two modulations were substantially different. We suspect 
that the difference came from the fact that the first-order approximation for the ellipsoidal effect is 
not accurate enough. 

The beaming effect increases (decreases) the brightness of any light source approaching (re- 
ceding from) the observer, e nabling the derivation of the stellar Doppler radial- velocity modulation 
from its precise photometry ( Loeb & Gaud 12003 : Zucker. Mazeh & Alexander"ll2007 ). The ampli- 
tude of the Doppler beaming is on the order of 4V re i/c, where V re i is the radial veloc ity of the source 
relative to t he observer a nd c is the speed of light (e.g.. Rvbicki & Lightmanll979 ). The ellipsoidal 
effect (e.g.. lMorriJl985L due to the tidal distortion of the star by the gravity of its companion, is a 
well known effect in short-period binaries. 

Given the stellar mass and radius and the or bital period of KOI- 13, the expected amplitudes 
of the beaming and the ellipsoidal effects (e.g., Fa igler & Mazehll201 II) are ~l.5(m p /Mj llp ) and 
~l5(m p /Mj up ) ppm. If indeed the secondary of KOI-13 is a brown dwarf, then the beaming ampli- 
tude is expected to be on the order of 30 ppm or even larger and the ellipsoidal effect should be ten 
times larger. Such periodic modulations, although much smaller than in the aforementioned case 
of KIC 10657664 , can be detected in the Kepler data, the ellipsoidal effect in particular. 

The two effects were already detected in a similar case — CoRoT-3, where a bro wn dwarf with 



a mas s of 22 Mj up was discovered to orbit an F-type star with a period of 4.3 days dDeleuil et al 



2008). In that case 



Mazeh & Faigler (2010) detected the beaming and the ellipsoidal effects with 



amplitudes of 27 + 9 and 59 ± 9 ppm, respectively, consistent with the mass previously known 
from the RV measurements. There is no reason why the Kepler publicly available data, which 
has a similar time span to that of CoR oT-3, would not reveal periodic modulations with similar 
amplitudes. In fact. IWelsh et a l. (2010) identified in the Kepler data the ellipsoidal effect induced 



by the known planet HAT-P-7b. 



Szabo etal 



(201 1) presented in their analysis of KOTO another smooth periodic modulation, 
well observed in short-period binaries — the reflection/emission modulation (referred to hereafter 
as the reflection modulation), induced b y the stellar light reflec ted by the day side of the secondary 
toge ther with its the rmal emission (e.g. Jcowan & Agollboi ll . In the folded light curve presented 
bv lSzabo et alJ (1201 ll) the occultation of the secondary by the primary star was also a very promi- 



nent feature. Given the secondary radius, the reflection phase modulation and the occult ation depth 
can be used to put some constrains on the secondary albedo and temperature (e.g., 



Snellen et al 



2009; 



Demorv et al 



201 11) . The derivation of these parameters is another opportunity to explore 



the nature of KOI-13.01. 

In this paper we present our analysis of the Q2 publicly available Kepler data of KOI- 13, detect- 
ing and deriving the amplitudes of the beaming, ellipsoidal and reflection modulations (Section 2). 
After this paper was submitted, Q3 data were released, so we repeated the analysis with the Q3 
data, the results of which are presented in Section 3. The results of the two quarters are quite simi- 
lar. We use the amplitude of the beaming modulation to estimate the mass of KOI-13.01 and show 
that it is probably a massive planet (Section 4). We use the derived amplitude of the reflection and 
the occultation depth to discuss the albedo and temperature of KOI-13.01 (Section 5). Section 6 
summarizes our findings. 
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2. Data analysis — Q2 data 

The analysis we present in this section is based on the long cadence, raw Q2 Kepler data, spanning 
88.9 days, from June 2009 till September 2009. These data were publicly available when we sub- 
mitted the paper. We decided not to use the Q0 and Ql data, as the transits of KOTO in that part 
of the light curve seemed inconsistent with the rest of the publicly available data. 



2.1. Preprocessing the light curve 



We 'cleaned' the data in the same way iMazeh & Fai gler (2010) preprocessed the CoRoT-3 light 
curve. This includes: 



- Transits removal. Measurements taken during and around the 50 transits of KOTO were re- 
moved from the analysis (but see below the analysis of the transits themselves.) 

- Dilution correction. We su btracted from all d ata points 0.45 times their median, which is the 
dilution factor according to J 



Szabqetal 



(120111) . 

- Outlier removal. We identified 58 outliers by calculating the running median and r.m.s. of 31 
measurements around each point, and rejecting measurements that were 3<x (4cr) or more higher 
(lower) than their corresponding median. This left a total of 3555 data points. 



- Long -term detrending with a cosine filter. We used a discrete cosine transform ( Ahmed et al 



1974), adapted to the unevenly spaced data we had in hand (for details see IMazeh & Faiglei 



(2010)). We fitted the data with a linear combination of the first 51 low-frequency cosine func- 
tions, the last of which is with a period of 4.351 d, and subtracted those functions from the 
data. 



We took special measures to ensure that the results did not depend on the long-term detrending. 
First, we decreased the number of low-frequency cosine functions till all the derived amplitudes 
converged. This left us with only 51 frequencies removed by the detrending. Second, we tried 
an independent smoothing algorithm, with a time span of 3 orbital periods, which yielded almost 
identical amplitudes. We therefore concluded that our detrending is a robust step of the analysis. 

Figure 1 shows the light curve before and after the removal of the long-term trend. The r.m.s. 
of the cleaned light curve is 80 ppm. 



2.2. Fitting amplitudes for the beaming, ellipsoidal, and reflection effects 



We fitted the cleaned data with a model that included the ellips oidal, beaming and reflection effects 
(hereafter the BEER model, following Faigler & Mazehfeoi 1 ). We approximated the beaming and 
the ellip soidal modula t ions b y pure sine/cosine functions, using mid-transit, t tmn = 2454964.93 
HJD, of iBorucki et al.1 d201 II) as phase zero. The beaming effect was presented by a sine function 
with the orbital period, and the ellipsoidal effect by a cosi ne function with ha lf the orbital period. 
The reflection was approximated by the Lambert law (e.g., Demorv et alJboi 1 ). 



In this approximation we expressed the relative stellar flux modulation AF/ F, where F was the 
averaged detrended, corrected for dilution, stellar flux, as a function of the orbital phase 



2n , 
- ~ tlran > ' 



(1) 



4 



Mazeh et al.: KOI-13.01 — a massive planet 



x 10 



S 



£-2 




mm 



,x 1 0"' 



10 



20 



30 



40 



50 



60 



70 



80 



90 




10 



20 



30 



40 50 
HJD - 2455000 



60 



70 



80 



90 



Fig. 1. The Q2 raw light curve of KOI- 13, before (lower panel) and after (upper panel) detrending. 
The plotted flux is relative to the mean flux, corrected for the contribution of KOT13B. The long- 
term model is presented by the red line. Note the different scales of the two panels. 



where P or b was the orbital period. We then wrote 

AFbeam(0 



F 

AF e iiip(f) 

F 

AF refl (f) 



= A beam sin cf> , 

A e uip COS 2(f) , 

sin (p + (n - 4>) cos <p 



= -A 



re 11 



(2) 
(3) 
(4) 



where the coefficients, A e iiip, Ab eam ing and A re fl were all positive. Note that we defined A re fl to rep- 
resent half the peak-to-peak variation of the reflection effect, to be consistent with the other two 
amplitudes. 

We fitted the cleaned, detrended light curve of KOTO with a 5-parameter linear model, 

sin (p + (n - (p) cos cp 



A1 B eer(0 = ao + Abeam sin (f> - A eBip cos 2cf> - A refl 



+ ci2 S sin 2(f> . 



(5) 



This was similar to what has been done bv lMazeh & Faigleri (120101) and lFaigler & Mazehl (1201 lb . 
except for the expression of the reflection effect. It was important to include the first overtone sine 
function in the model as a check of our approach. We expected this term to be small, as it had no 
specific role in the model. The fitting process could find any value, positive or negative for the five 
parameters. However, we did expect Ab eam , A e ii; p and A re fl to be positive and «2s to be close to zero. 

Figure [2] shows the best-fit model and Table Q] lists the resulting amplitudes, after applying 
the dilution correction due to KOI-13B. The table includes also the results of fitting the model to 
Q3, as discussed in the next section. All three amplitudes come out to be positive, as expected. In 
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Table 1. The best-fit coefficients of the three effects and the occultation, using the Q2 and the Q3 



data 

Q2 Q3 

BEER model: 

A bcilra 8.6+1.1 10.4 ±1.0 ppm BEaming 

A e iiip 66.8 ± 1.6 69.1 + 1.5 ppm Ellipsoidal 

A refl 72.0+ 1.5 72.0 ±1.4 ppm Reflection 

a 2s -1.0 ±1.0 -3.3 ±1.1 ppm — 
Primary transit: 

a/R, 3. 16 ± 0.08 3. 17 ± 0.06 Semi-major axis to stellar radius 

r,,/R, 0.0908 ± 0.0005 0.0907 ± 0.0005 Planetary to stellar radii 

b 0.75 ± 0.01 0.75 ± 0.01 Impact parameter 

Secondary occultation: 

depth 163.8 ± 3.8 157.8 ± 3.5 ppm 

A</> -0.0007 ± 0.0005 -0.0012 ± 0.0005 



addition, the fourth coefficient is substantially smaller than the smallest amplitude, and is consistent 
with zero, as expected. The errors of the four coefficients are around 1-2 ppm, attesting to the high 
precision of the Kepler satellite. 

The ellipsoidal modulation, with two peaks, one at a phase of 0.25 and the other at 0.75, is 
very prominent in the figure. So is the reflection effect, which causes the folded light curve to be 
brighter at phase 0.5 than at phase 0. The beaming modulation, which causes the peak at 0.75 to be 
lower than the one at 0.25, is quite small but nevertheless noticeable even by eye, with statistical 
significance of almost ~8cr. 

One can also easily notice the secondary occultation at phase 0.5. We fitted the occultation 
with a plain trapezoid, using only two parameters — the depth and the phase of the occultation. It 
turned out that the shape of the trapezoid was determined by the three geometrical parameters of 
the systems — a/R*, the ratio of the semi-major axis to the stellar radius, R* / r p the ratio of the stel- 
lar to the planetary radii, and b, the impact parameter. Those parameters could have been derived 
accurately from the planetary transit. Therefore, although it was not the goal of this study, w e nev- 



ertheless fitted the long-cadence data of the planetary transit with the lMandel & Agoll d2002l) model 



implemented in MATLABU, using quadrati c limb darkening , with parameters that were found to 



best fit the data. Our fit used the value of the 



Mandel & Agoll (120021) model at the Kepler timing of 



each point, and ignored the finite sampling interval of the data. 

The best derived values of the three geometrical parameters, used to fit the secondary occulta- 
tion, are also listed in Table 1. The phase of the occultation is very close to <p = 0.5. The difference, 
denoted in the table as A0, is consistent with zero, probably indicating a highly circular orbit. 

2.3. Another unexpected periodicity 



In ord er to look for additional periodicities we calculated the Lomb-Scargle periodogram (IScargle 



1982) of the residuals, derived from the preprocessed data by removal of the BEER and the sec- 



ondary model. The residuals and their periodogram are plotted in Figure 3. 



Can be found now at www.astro.washington.edu/users/agol/transit.html 
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300- 




Fig. 2. The folded cleaned Q2 light curve of KOI- 13 with the best-fit BEER model. The residuals 
are plotted at the bottom of the figure, with a trapezoid fitting of the secondary occultation. 
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Fig. 3. The Q2 residuals of KOTO, after removing the BEER and the secondary model (upper 
panel) and their Lomb-Scargle periodogram (lower panel). 
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Fig. 4. The Q2 residuals of KOI- 13, after removing the BEER and the secondary model, folded with 
the newly found period of 1.0595 d, and binned into 200 bins. The error bar of each bin presents 
the scatter of the points in that bin, divided by the square root of the number of points in that bin. 



The figure shows a prominent peak at a frequency of 0.9438 + 0.001 d _1 , corresponding to a 
period of 1 .0595 d. Although suspiciously close to the notorious 1 -day period, which might point to 
some telluric origin, we can not see any reason to exclude an astrophysical modulation (see further 
discussion in the next section.) The folded data with the newly found period is shown in Figure 4, 
which shows that the amplitude of the modulation (peak-to-peak) is ~60 ppm. 



3. Analysis of Q3 data 

After we submitted this paper, Kepler Q3 data were released, so we repeated the analysis with the 
newly available raw data. We opted to analyze the two data sets separately, in order to have an 
independent handle on the errors of the derived parameters. 

The Q3 light curve includes data taken during 89.3 days, from September 2009 to December 
2009. The preprocessing analysis was similar to the one of Q2, after which we were left with 3613 
data points with an r.m.s. of 77 ppm. We then fitted the BEER and the trapezoid models, the results 
of which are listed in a separate column of Table 1 . 

The newly derived parameters of the three modulations are quite similar to the ones of Q2 data. 
All differences are within 1.2cr of their corresponding errors. Even the difference in the secondary 
occultation depth, 6.0 + 5.3 ppm, is not significant. Therefore, for the following discussion we 
adopted the values derived from Q2. 

In order to check the additional periodicity of 1 .0595 d found in Q2, we plotted the power- 
spectrum of the residuals found in the Q3 data. The derived spectrum shows a few more peaks, 
as can be seen in the upper panel of Figure [5] The fact that those additional peaks were not seen 
in Figure |4] might attest to the better quality of the Q3 dataset. The additional frequencies can be 
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Fig. 5. Upper panel: The power spectrum of the Q3 residuals of KOI- 13, after removing the BEER 
and the secondary model. Lower panel: The same, after removing the main periodicity of 1 .0595 d 
and its four overtones. For comparison, we kept the same scale for the two panels. 

seen better in the lower panel of the figure, which shows the Q3 power spectrum after removing 
the 1.0595 periodicity with its four overtones. 



4. The mass of KOI-13.01 

In this section we derive the mass of the planet candidate KOI-13.01, based on 

- the amplitude of the beaming and ellipsoidal effects, 

- th e first order theoretical e xpectatio n of the beaming and el lipsoidal effects, as parameterized 
bv lMazeh & Faigleri (12010 ) (see also 

- the stellar parameters, as derived by 



Faigler & Mazeh 



Szabo et al 



(2011). 



2011), and 



Using Faigler & Mazehld201 II) approximations and notations we can derive the secondary mass 



from the beaming and ellipsoidal amplitudes: 



^2,beam 



«2,ellip - 



0-37 A beam / M* 
«beam lppm\M 



2/3 



1 day 



1/3 



M 



3 up 



0-08 Aellip [R*\ (M,y ( forb V 

ldayj 



a e iiip lppm \R, 



Mr. 



M 



J up 



(6) 



(7) 



In the above expressions M„ and R t are the primary mass and radius, respectively. The two ni2's 
present the estimated mass of the secondary from the two amplitudes, and the two or's are numerical 
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Table 2. KOI- 13 A parameters, as derived by J 
Parameter Derived value Unit 



Szabo etal 



(2011) 



M, 
R, 
T, 



2.05 M Q 
2.55 R B 
8500 ±400 K 



Dilution factor 1.82 



factors which are on the order of unity. For completeness, we list the relevant parameters of the 
system in Table 2. The period value is 1.763582 + 0.0000014 d. 

The abeam factor includes one component that ori ginates from the fact that the stellar spectrum 
is Doppler shifted. To estimate this factor we followed 



shifted spectra from the l ibrary of 



Castelli & Kurucz 



Faigler, M azeh et al.l(l201 II) and numerically 



(2004) models, which were then integrated 



over the Kepler bandpass ( Zuck er. Mazeh & Alexander 1120071) . yielding for each spectrum an a\, e , 
factor. The value we adopted for KOI- 13, 



«be. 



0.63 + 0.06 



(8) 



was derived by interpolation bet ween the available models of the library that were close to the 
temperature of KOI- 1 3 given by ; 



Szabo et al 



d201 ll) . We estimated the uncertainty by calculating 
the interpolated abeam values within the T e ff error range, and by taking into account the possible 
gravity and metallicity ranges. All together, the error on abeam amounted to 10%. 



The uncertainty of g e iiip was much larger. We ran numerous models of EBOP ( Etzel 



1980; 



Popp er and Etzell 1 1 98 lb . EBAS dTamuz et alj 120061) . and Wilson-Devinney dWilson & Devinnev 
197 lh codes to fit the KOI- 13 lightcurve and tried to estimate the value of a e ii; p and get a better 
value than the unity default estimate, to no avail. This might be due to the high gravity darkening 
and the fast rotation of KOI- 13. We therefore adopted a conservative estimate of the numerical 
factor of 



*ellip 



= 1 



+ i 



(9) 



To complete the error est imate of the s econd ary mass we needed the uncertainties of the stellar 
mass and radius of KOI- 13. 



Szabo et al 



(2011) did not give error estimates for their values and 



therefore we estimated the stellar mass and radius error to be 10%. 
Using the values of Table 2, we got: 



W 2 ,beam -10 + 2 M Jup m 2; ellip - 4^ M. 



J up 



(10) 



Due to its large fractional error, the mass derived from the ellipsoidal modulation was not very 
useful. We therefore adopted the value derived from the beaming modulation. This value indicated 
that the secondary mass was very probably in the higher mass range of the planetary regime. 



5. The albedo and temperature of KOI-13.01 

The strong insolation of the A-type star that hits KOI-13.01 should heat up its day side to a rela- 
tively high temperature. If one assumes that the planet has reached a full synchronization, that all 
stellar energy is absorbed and no transfer of energy to the night side takes place, then the day-side 
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Fig. 6. The difference between the occultation depth and the phase modulation, Ans, as a function 
of a/R*, based on Q2 and Q3 data, combined. The vertical lines present the best-fit a/R* and its lcr 
error estimation, derived from the planetary transit. 

temperature can be of the order of 3500 K. This is high enough to explain the observed phase mod- 
ulation with peak-to-peak amplitude of 144 + 1 ppm. Obviously, according to this approach the 
night side does not emit any light and therefore the occultation depth should be equal to the phase 
modulation amplitude. The depth values we derived in the previous sections — 163.8 + 3.8 for Q2 
and 157.8+3.5 for Q3, indicate that the occultation depth might be slightly deeper. According to the 
simplistic approach that assumes no light from the night side, the small derived difference between 
the occultation depth and the phase modulation, denoted here as A^s > is necessarily insignificant. 

However, if the observed depth difference is real, we have to assume that it comes from the 
night side, and therefore it gives us a rare opportunity to explore the night-side thermal emissivity, 
the Bond and geometrical albedos of the planet, and finally the effective temperatures of the day 
and the night sides of the planet. 

To assess the significance of the non-zero detection of A^s we first note that the figures of Q2 
and Q3 in Table 1, taken at face value, yield A^s = 16.8 + 2.8 ppm. However, as mentioned in the 
previous sections, these values depend on the geometrical parameters of the system, derived here 
by the analysis of the transit. In particular, the trapezoid depth depends on a/R*, the ratio of the 
semi-major axis to the stellar radius. To further estimate the reality of the A^s detection we plot in 
Figure|6]the weighted average value of A^s, obtained from Q2 and Q3, as a function of a/R*. The 
figure also shows the lcr error on a/R*. From the figure we conclude that our best estimate for the 
night-side emissivity is 

A NS = 16.8 +4.5 ppm. (11) 

The figure suggests that the detection of the small difference between the occultation depth and the 
phase modulation might be real. 

We now move to exploit the value derived for A^s for better understanding of KOT 13.01 energy 
budget. To begin, we find that the night-side effective temperature, 7n, needed to account for the 
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observed A NS is 2600+150 K, assuming r p /R, = 0.0908 + 0.0003 (see Table 1). In the error budget 
we took into account the uncertainties of the stellar temperature, the stellar and planetary radius 
ratio, and the error of A^s ■ The derivatio n of assumes a blac k-body emissivity of the night side, 
obviously a very simplistic assumpti on dCowan & Agolll201 lh . A detailed model of the planetary 
atmosphere (e.g., Fortne v et al is well beyond the scope o f this paper. 

To proceed, we use the formalism of ICowan & Agoll (1201 ll) . They approximate the day-side 
effective temperature as 

(RA K2 /„ c \i/4 



and that of the night side as 



(12) 



(13) 



In these expressions T e g jt is the stellar effective temperature. Two parameters describe the energy 
budget of the planetary day side — Ab is the "effective" Bond albedo, which represents the fraction 
of the total energy of the star reaching the planetary surface that is not re-emitted thermally by the 
planetary atmosphere, and e is the fraction of the energy absorbed by the day side and re-emitted 
thermally that is transferred to the night side. The estimation of the night-side temperature gives 
us a strong constraint on the day-side temperature and the fraction of the energy transferred to the 
night side. This is done by writing 

/8 5\ 1/4 

7W N (---) , (14) 

and using the derived to calculate the dependence of 7b on e. We can then derive the Bond 
albedo for a given To- This is depicted in Figure [7] We also plotted in this figure t he geo met- 



rical albedo, A„, as a function of I'd, as was done, for example, by 



Demoryetal 



Santerne et al. 



(2011) and 



d2011l) . and the corresponding e value. 
The figure sugges ts that if indeed the de rived 7n is correct, Ab is substantially larger than A g in 
the Kepler band (e.g., Cowan & Agollboi 1 ) for most cases. 



6. Discussion 



The analysis we present here, based on the beaming modulation detected in the Q2 and Q3 publicly 
available data, suggests that KOTO. 01 is a massive planet, with a mass in the range of 8-12 Mj up . 
The mass derived from the ellipsoidal modul ation is not very useful, as the expect ed amplitude is 
derived from first-order approximations (e.g. 



Morris 1985 



an 



Morris & Naftilan 1993), and therefore 



2011) is not 



the degree of its applicability to hot rapidly rotating stars like KOI-13A (Szabo et al 
clear. When more stars with beaming and ellipsoidal modulations are detected ( Carter et al.ll2.Ql ll) 
we expect the theory to be better tuned and become more accurate. 

As of July 201 1, KOI- 13.01 has the largest radius known for an exoplanet. To locate its position 
in the radius-mass parameter plane we plot the known transiting planets in Figure [8] showing that 
KOTO. 01 has indeed the largest known radius. It seems as if the upper envelope of the planetary 
radius reaches its maximum slightly below Jupiter mass. If our estimate for its mass is correct, 
KOTO. 01 is substantially larger than planets with similar masses, but its density is similar to that 
of Jupiter itself. Much theoretical work has been done to explain the large radii of some planets 



12 



Mazeh et al.: KOI-13.01 — a massive planet 




2600 



2800 



3000 



3200 



3400 



3600 



day 



Fig. 7. The effective Bond albedo, Ag, the geometrical albedo, A g , and e, the energy fraction trans- 
ferred to the night side, as a function of the day-side temperature. One-sigma errors are plotted 
around each curve. 



(e.g., 


Guillo 


t & Showmanl2002 


Baraffe et alj|2004 


Chabrier & Baraffel2007 


Fortnev et al. 


2007; 


Miller et al. 


2009 


). It would be interesting to find observational indications whether KOI-13.01's 



Enoch et al 



20 lib and/or the 



Szabo et al 



large radius is associated with its large if not largest insolation (e.g., 
relativ ely large mass of its host star and/or the young age of the system, estimated by 
J201 ll) to be 0.7 Gyrs. 

We have discovered additional periodicities in the light curve of KOI- 13, with amplitudes of 
30 ppm or less and frequencies in the range of 0.5-2 d _1 . We suspect that the origin of these 
modulations is stellar. 

The analysis suggests that the secondary occultation depth might be deeper than the reflection 
peak-to-peak modulation. If tr ue, this is a way to es timate the thermal emission of the planetary 
night side of KOTO. 01 (e.g.. ICowan & Agolll201 lb . We estimate its effective temperature to be 
2600 ± 150 K, assuming the detected difference is real, and using a simplistic black-body emis- 
sivity approximation. We derive the planetary geometrical and effective Bond albedo, the latter 
representing the fraction of the total energy of the star reaching the planetary surface that is not 
re-emitted thermally by the planetary atmosphere, as a function of the day-side temperature. Our 
analysis suggests that the Bond albedo is probably substantially larger than the geometrical albedo. 
This might indicate that some of the stellar photons not seen in the Kepler band do not share the 
same processes as the optical photons. For example, UV photons, beyond the Balmer limit, might 
be absorbed by Hydrogen atoms at the first excited state in the upper layers of the atmosphere and 
their energy might not be directed to heat the planetary atmosphere. 

Unfortunately, the detection of the night-side thermal emission was not highly significant. With 
Kepler additional data, beyond Q2 and Q3, one should be able to determine with higher significance 
and precision the magnitude of the beaming, ellipsoidal and reflection effects, and compare the 
secondary depth with the phase modulation amplitude. 
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Fig. 8. The radii of the transiting planets as a function of their masses, taken from 
http://www.exoplanet.eu KOI-13.01 is plotted as a red square. The dashed lines are constant den- 
sity lines, of 0.5, 1 and 2 Jupiter density. 



KOTO is one of the very few A-type stars known to host a hot Jupiter (IDesort et al.ll2010h . 
More cases like KOTO might present a unique opportunity to shed some light on the popula- 
tion of giant planets around stars more massive than most of the known planet-hosting stars, and 
compare them to the planets around G- and K-type stars. Right now this is done mainly by RV 



Bowler et al 



search of a sample of sub-giants that are believed to be evolved A-type stars (e.g., 
2010). Those sub-giants might not be able to retain hot Jupiters with short periods like that of KOI- 
13.01. We assume that many non-eclipsing KOI- 13. 01 -type planets are lurking in the Kepler and 
CoRoT data. In fact we are searching in both data s ets for such systems with the BEER algorithm 
( Faigler & Mazehlboi 1 : Faieler. Mazeh et alJEoil ). 

The estimated mass of KOI-13.01 might justify a recognition of its planetary nature, even with- 
out c onfirmation by fo llow-up RV observations. This w as done already, for example, for Kep ler- 



9d (ITorres et al 



201 II) . Kepler- 10c (IFressin et al.ll2.01 II) . and Kepler- llg (iLissauer et al 



2011b . It 



seems as if the recent large number of transit detections, especially of low-mass planets and early- 
type host stars, pushes us towards a new terminology of confirmed planets. 
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